Effects of initial state fluctuations on jet energy loss 
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The effect of initial state fluctuations on jet energy loss in relativistic heavy-ion collisions is studied 
in a 2+1 dimension ideal hydrodynamic model. Within the next-to-leading order perturbative QCD 
description of hard scatterings, we find that a jet loses slightly more energy in the expanding 
quark-gluon plasma if the latter is described by the hydrodynamic evolution with fluctuating initial 
conditions compared to the case with smooth initial conditions. A detailed analysis indicates that 
this is mainly due to the positive correlation between the fluctuation in the production probability 
of parton jets from initial nucleon-nucleon hard collisions and the fluctuation in the medium density 
along the path traversed by the jet. This effect is larger in non-central than in central relativistic 
heavy ion collisions and also for jet energy loss that has a linear than a quadratic dependence on its 
path length in the medium. 

PACS numbers: 12.38.Mh, 24.85. +p; 25.75.-q 



J2' 

9-; 



> 
O 
oo 

OV 
(N 

OO 

o 

(N 



X 



I. INTRODUCTION 



The observation of jet quenching through the suppres- 
sion of large transverse momentum single hadron, di- 
hadron, and 7-hadron spectra in relativistic heavy ion 
collisions [H, Q is one of the most important evidence for 
the formation of a strongly coupled quark-gluon plasma 
(QGP) in these collisions. Jet quenching is a measure 
of the energy loss of an initial leading jet parton as it 
traverses through the produced dense matter via mul- 
tiple scatterings Q. Theoretical studies on parton jet 
energy loss have concentrated on both gluon radiation 
induced by multiple scattering and elastic collision en- 
ergy loss. Due to the non-Abelian Landau-Pomcranchuk- 
Migdal (LPM) interference effect the radiative en- 
ergy loss shows a quadratic path-length dependence , 
which is in contrast to the linear path-length dependence 
of the elastic collision energy loss ficil - [l2j . Also, a cubic 
path-length dependence of the jet energy loss has been 
found in the strongly coupled limit of the QCD medium 
using the AdS/CFT correspondence [il^flilj. 

The study of jet quenching in heavy ion collisions has 
been carried out in the 1+1 dimension Bjorken hydro- 
dynamics (T§4l8| as well as the 2+1 and 3+1 dimension 
ideal and viscous hydrodynamics [l9l - l22l |. In these stud- 
ies, the initial conditions for the hydrodynamical evolu- 
tion were taken to be smooth in space. Recently, the ef- 
fect of initial event- by- event fluctuations on jet quenching 
has been investigated in the 1+1 Bjorken hydrodynam- 
ics [23| . It was found that the strong correlation between 
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the fluctuation in the spatial distribution of initial hard 
scatterings from which jets are produced and the fluc- 
tuation in the density distribution of the initial medium 
has significant effects on jet quenching. In particular, 
the jet energy loss is reduced after the inclusion of ini- 
tial fluctuations. However, the transverse expansion of 
the produced hot dense medium has been neglected in 
this study. In the present paper, we include the trans- 
verse expansion in studying the effect of initial fluctua- 
tions on jet energy loss by using the 2+1 dimension ideal 
hydrodynamic model of Refs. @, H^. For calculating 
the hadron spectra at large transverse momentum, we 
use the next-to-lcading order (NLO) perturbative QCD. 
Our results show that including the transverse expansion 
of the medium slightly enhances the energy loss of jets, 
contrary to the reduced jet energy loss found in Ref. [23| 
without the transverse expansion. We further investigate 
the effect of initial fluctuations for different path-length 
dependence of jet energy loss in the medium. 

This paper is organized as follows. We first give a brief 
description of the the 2+1 dimension ideal hydrodynamic 
model in Sec. Inland the jet quenching models in Sec. IIIII 
Results from our study are shown in Sec. IIVI We then 
present some discussions in Sec. |V]and finally summarize 
our study in Sec. [Vll 



II. 2+1 DIMENSION HYDRODYNAMICS 

In the 2+1 dimension ideal hydrodynamics, which as- 
sumes the boost invariance along the longitudinal direc- 
tion, the energy-momentum tensor T^ u and pressure p of 
a system can be expressed in terms of the proper time r 
and the two transverse coordinates x and y perpendicular 
to the beam direction p3 - l27j . Conservations of energy 
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and momentum then give 

d T (TT 00 ) + d x (TT 0x ) 
d T ( T T 0x ) + d x (rT xx ) 



-d y (rT°y) = -p, 
d y (rT xy ) = 0, 



d T {rT 0y ) + d x {rT xy ) + d y {TT yy ) 



0. 



(1) 



To solve these equations requires information on the ini- 
tial conditions of a collision, particularly the initial en- 
tropy density, and the equation of state of the produced 
matter. For the initial entropy density, it is taken as 



drj 



(i-0 



^part 
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where n par t and n co n are the number densities of partic- 
ipants and binary collisions, respectively. 

In heavy ion collisions, the initial conditions vary from 
event to event as the positions of colliding nucleons are 
randomly distributed according to the density distribu- 
tions of the colliding nuclei. Two nucleons are consid- 
ered as participants and a binary collision takes place at 
their middle point if the transverse distance between a 
nucleon from one nucleus and a nucleon from the other 
nucleus is shorter than <Ji n /ir, where a ln = 42 mb is the 
nucleon-nucleon inelastic cross section at RHIC energies. 
A smearing parameter a is then introduced in evaluating 
the number densities of participants and binary collisions, 



^part(coll) ( r ) — 
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H exp 
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(3) 



where is the transversal position of a participant (bi- 
nary collision). Here we use the same smearing parame- 
ter a for both the participant and binary collision num- 
ber densities. In the present study, we consider the two 
cases of a = 0.4 fm and 0.8 fm. Also, we choose the 
initial thermalization time to = 0.6 fm/c for starting the 
hydrodynamical evolution. 

For the equations of state, we use the quasi-particle 
model based on the lattice QCD data for the QGP and 
the resonance gas model for the hadron gas as in Refs. 0, 
l25l [29j . This model thus assumes the presence of a first- 
order phase transition and the critical temperature T c is 
170 MeV. We solve the hydrodynamic equations Eq. (TTJ) 
numerically by using the HLLE algorithm [30l - l32l | . 

The parameters C and £ in Eq. ([2J are determined from 
fitting the centrality dependence of the final charged- 
particle multiplicity (28|. Using the Cooper- Frye freeze- 
out formula and assuming that the multiplicity does not 
change after chemical freeze-out at temperature T = 160 
MeV 0, we obtain C = 19.3 and £ = 0.11. 

In studies with smooth initial conditions, both the par- 
ticipant number and the binary collision number densities 
are obtained from the thickness functions of the collid- 
ing nuclei evaluated from their density distributions. In 
the present study, they are obtained by averaging over 
a large number of initial fluctuating events. Because of 



the smearing parameter introduced in generating the ini- 
tial conditions for hydrodynamical evolutions, the result- 
ing smooth participant number and the binary collision 
number densities have a larger spread in space than that 
obtained from the nuclear thickness functions. 



III. JET QUECHING MODELS 

For a jet of energy E produced at the position r from 
a hard nucleon-nucleon collision and moving along an 
azimuthal angle 4> in the transverse plane of a nucleus- 
nucleus collision, its total energy loss can be expressed 
as 



AE = / dTf{E,t/,,T,T)p(T,4>,T), 



(4) 



where p(r, r) is the local parton density at time r along 
the jet path, and the function f(E, <j>, r, t) is the jet en- 
ergy loss per unit time through a unit density of medium. 

Averaging over the creation positions and moving di- 
rections of the jet in the transverse plane gives 



(AE) = i- J d^d 2 rdTn(r)f(E l( f>,r,T)p(r^ 



(5) 



where n (r) = n co ii(r)/iVbin, with n co n(r) and iV co n denot- 
ing, respectively, the number density of binary collisions 
at r and the total number of binary collisions, is the prob- 
ability density for jet production at r. The average jet 
energy loss rate along the jet path is then 



J d(t>d 2 m(r)f(E,4>,r,T)p(r, l 



(6) 



d(AE) 
dr 

According to recent theoretical studies [15|, UjJ, I33H3 
the total quark energy loss in a finite and expanding 
medium is approximately given by 

'dE\ r dT f^ro\ a p(r,r) 
dL 



AE 



TO 



p(r ,0) po 



(7) 



In the above, a is a parameter with possible values 0, 1 
and 2, corresponding, respectively, to linear, quadratic, 
and cubic path-length dependence for the jet energy loss; 
and p^ and are, respectively, the four momentum of 
the jet and the four flow velocity of the local medium. 
The average energy loss per unit length (dE/dL) has the 
following parametrization [34|: 



dE\ 

— ) = eo(E/p 
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(8) 



where eo is the energy loss parameter with a value that 
is 9/4 times larger for a gluon than for a quark, and p is 
the Dcbye mass. Their values arc determined from fitting 
the experimental data for the nuclear modification factor 
in the most central A + A collisions using the smooth 
initial conditions. This leads to the following jet energy 
loss rate through a unit density of medium: 



/(£,</>, r,T) = 



dE 
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IV. RESULTS 



DISCUSSIONS 



In this Section, we present results for the nuclear mod- 
ification factor of jets in both central and mid-central 
Au+Au collisions at ^/snn = 200 GeV. For a given event 
i, the nuclear modification factor is defined as 



R 



AA 



da AA /dp^dy 



(10) 



where the hadron spectra in heavy-ion collisions 
du % AA j dp\dy are calculated from the NLO pQCD with 
modified fragmentation functions due to jet quenching 
flol Il6j . Specifically, the cross sections for the hard scat- 
tering are calculated using the CTEQ6M parameteriza- 
tions [36| for the parton distributions in a nucleon and 
including both 2 — > 3 tree level contributions and 1-loop 
virtual corrections to 2 — > 2 tree processes [13] ■ Further- 
more, the AKK08 parameterizations [38| are used for the 
parton fragmentation into hadrons. 

The nuclear modification factor for the case of the 
fluctuating initial conditions is then given by the aver- 
age of Eq. (fTU]) over all events. For the case of smooth 
initial conditions, the nuclear modification factor is ob- 
tained, on the other hand, from the ratio of correspond- 
ing charged hadron spectrum to that of the proton- 
proton collision multiplied by the average number iV co ii 
of nuclcon-nuclcon collisions. 
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FIG. 1: (Color online). Nuclear modification factors with and 
without initial fluctuations in 0-5% and 40-50% centralities 
of Au+Au collisions at y/s = 200 GeV, respectively. The 
experimental data are taken from [39| . 



Fig.[T]shows the nuclear modification factors of high-pr 
particles in 0-5% and 40-50% centralities of Au+Au col- 
lisions at y/s = 200 GeV with and without initial fluctu- 
ations. It is seen that including initial fluctuations leads 
to a smaller Raa, an d the effect is stronger for larger 
initial fluctuations (corresponding to smaller a) and in 
noncentral than in central collisions. 



The smaller Raa in the case of fluctuating initial con- 
ditions obtained in the present study is opposite to the 
result reported in a previous study based on the 1+1 
boost invariant hydrodynamics [23[, where a larger Raa 
was obtained when initial fluctuations were included. To 
understand this difference, we define the jet energy loss 
difference 6(AE) = (AE) auc - (AE) smth between the av- 
erage energy loss calculated with fluctuating initial condi- 
tions and that with smooth initial conditions. From Eq. 
©, we then obtain the following rate for this difference 
along the jet path: 



d[6(AE)} 



1 

2^ 



d(f)d 2 rf{E,(j),r,T){dnSp), (11) 



where Sn and Sp are, respectively, the differences in the 
jet production probabilities and the medium densities in 
the cases of fluctuating and smooth initial conditions. 
We note that the terms pSn and n5p are not present 
in Eq. (|lip since both vanish after integration over the 
jet production positions and moving directions. Eq. (|11[) 
indicates that the energy loss difference is determined not 
only by the fluctuation in the jet production probability 
density but also by the fluctuation in the local parton 
density on the jet path [23| . 

200GeV AuAu->ti° 0-5% ct = 1 p T =20GeV 
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FIG. 2: (Color online). Rate of jet energy loss difference, 
averaged over all jet paths, in 0-5% central Au + Au collisions 
at y/s — 200 GeV for smearing parameter a — 0.4 fm and 
0.8 fm. The inset shows the averaged energy loss rate along 
the jet path calculated with fluctuating and smooth initial 
fluctuations. 

In Fig. [2J we show the rate of the energy loss difference 
for jets with the transverse momentum px = 20 GeV, 
averaged over all jet paths, in 0-5% central Au+Au col- 
lisions at y/s = 200 GeV for smearing parameter a = 0.4 
fm and 0.8 fm. It shows that the correlation between 
the fluctuation in the production probability of initial 
parton jets and the fluctuation in the local medium den- 
sity is positive during the initial stage of jet propagation 
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but changes to negative in the later stage, resulting in 
enhanced and reduced energy losses, respectively. This 
result is consistent with that in Ref. [23| for a trans- 
versely static medium. The net effect of the fluctuations 
on the jet energy loss is determined by the sum of the 
positive and negative differences. As shown in the inset 
of Fig. [2J which gives the averaged energy loss rate along 
the jet path calculated with fluctuating and smooth ini- 
tial fluctuations, most energy losses happen close to the 
initial path of the jet. Because of the dominance of the 
initial positive difference, the total energy loss calculated 
with fluctuating initial conditions is greater than that 
with smooth initial conditions. 

The relation between the jet propagation and the 
medium evolution can be further clarified if we approx- 
imate the time evolution of the parton density p(r, r) 
along a jet path in the hydrodynamic evolution of the 
medium as p(r,r) ~ I/t* 9 . According to Eq. ©, the 
time evolution of the medium-dependent jet energy loss 
can then be simply written as 



AE 



r a-/3 



(12) 



if we neglect the small flow effect. As the jet transverses 
through the medium, its increasing energy loss with the 
path- length (r Q ) is thus suppressed by the decreasing 
density of the bulk medium (1/r' 3 ). Consequently, the 
total effect of the initial fluctuations on jet quenching is 
related to the competition between the path-length de- 
pendence and the medium-density evolution dependence 
of the jet energy loss. Since a is always smaller than /3 in 
our study, the total energy loss mainly takes place during 
early times when the correlation is positive, thus resulting 
in more energy loss in the case of fluctuating initial condi- 
tions. In non-central collisions, the fireball expands faster 
than in central collisions (^j 4 °- 50 % > (3°~ 5% ) y so most en- 
ergy loss happens earlier than in central collisions. As a 
result, the Raa for the case of initial fluctuations with 
a = 0.4 fm decreases by 4% in central collisions and by 
8% in the 40-50% centrality of the collisions as shown in 
Fig.ffl 

For the longitudinal expanding medium and the radia- 
tive energy loss mechanism considered in Ref. [23j . j3 is 
much smaller than a, so the total energy loss is domi- 
nated by the contribution during later times when the 
correlation has large negative values. Therefore, the en- 
ergy loss is smaller and the Raa obtained in Ref. [23| 
for the fluctuating initial conditions is larger than in the 
smooth case. 

The path- length dependence of jet energy loss depends 
on the energy loss mechanism. It is linear (a=0) for elas- 
tic energy loss, quadratic (a = 1) for radiative energy 
loss p-Qjll, and cubic (a = 2) for energy loss based on 
AdS/CFT for the strongly coupled QCD [H03. How 
the different power of path-length dependence affects the 
effect of initial fluctuations on jet energy loss is an inter- 
esting question. We illustrate this effect by considering 
that has either linear or quadratic or other path-length 
dependence but with the energy loss parameter eo in Eq. 



([HI) always determined from the smooth case by fitting 
the experimental data of the nuclear modification factor 
in the most central A+A collisions. 
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FIG. 3: (Color online). Rate of jet energy loss difference av- 
eraged over all jet paths in 0-5% central Au+Au collisions for 
different powers of the path-length dependence of jet energy 
loss, a = 0, 1 and 1.5. The inset is the average energy loss 
rate along the jet path for a — and 1 
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FIG. 4: (Color online). Effects of initial fluctuations on jet 
quenching manifested by the nuclear modification factor of 
Pt — 7 — 9 GeV hadrons as a function of the power index 
of the path-length dependence of the jet energy loss for three 
different definitions of the medium density in central Au+Au 
collisions. 

In Fig. [3J we show the rate of jet energy loss av- 
eraged over all jet paths in 0-5% central Au+Au colli- 
sions for different powers of the path-length dependence 
of jet energy loss, a = 0, 1 and 1.5. For the linear path- 
length dependence (a = 0) of jet energy loss, most en- 
ergy loss takes place in the initial positive correlation 
region of the jet path, so the initial positive correlation 
dominates the fluctuation effect. For the quadratic path 
length dependence (a = 1), the peak for the jet energy 
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loss rate is shifted closer to the negative correlation re- 
gion as shown in the inset of Fig. [31 so the dominance 
of initial positive correlation is weakened by the nega- 
tive correlation. Therefore, the energy loss for the linear 
path-length dependence is greater than the energy loss 
for the quadratic path-length dependence. This conclu- 
sion is supported by the results shown in Fig. 0] for the 
ratio (R A U £ - R S AA h )/ R AA h of p T = 7 - 9 GeV hadrons 
in central Au+Au collisions as a function of a. For the 
smaller smearing parameter a =0.4 fm, the initial fluc- 
tuating conditions decrease the suppression factor by 8% 
for the linear path-length dependence of jet energy loss 
while by 4% for the quadratic path-length dependence 
of jet energy loss. The fluctuation effect with the larger 
smearing parameter a = 0.8 fm is smaller than that with 
a = 0.4 fm. 



VI. SUMMARY 

Based on the 2+1 dimension ideal hydrodynamics, we 
have studied the effect of initial fluctuations on jet en- 
ergy loss in relativistic heavy-ion collisions within the 
description of the NLO pcrturbative QCD. Our results 
show that fluctuating initial conditions lead to slightly 
more energy loss than smooth initial conditions. In gen- 
eral, the jet energy loss increases with time due to its 
path-length dependence but this increase is suppressed 
by the decreasing medium density with time. Where the 
total energy loss mainly takes place along the jet path is 
determined by the competition between the path-length 
dependence of jet energy loss and the time dependence 
of the medium density. For fluctuating initial conditions, 
our results for the rate of the average energy loss dif- 
ference between the two cases of fluctuating and smooth 
initial conditions show that the correlation between the 



fluctuation in the production probability of initial par- 
ton jets and the fluctuation in the local medium density 
is positive during the early times along the jet path and 
negative during the later times. Consequently, the net ef- 
fect of initial fluctuations on jet energy loss is determined 
by whether the energy loss mainly takes place when this 
correlation is positive or negative. The total energy loss 
in the fluctuation conditions is then larger than that in 
the smooth case if most energy loss takes place when 
the correlation is positive, while it is smaller if it takes 
place when the correlation is negative. Our results fur- 
ther show that the initial positive correlation dominates 
the fluctuation effect for linear and quadratic path-length 
dependence of jet energy loss in central as well as in non- 
central A+A collisions. However, because this dominance 
is stronger in non-central collisions than in central col- 
lisions, the difference between the nuclear modification 
factors calculated with fluctuating initial conditions and 
smooth initial conditions in non-central A+A collisions 
is greater than that in central A+A collisions. Similarly, 
the jet energy loss for the linear path-length dependence 
is more affected by the fluctuation effect than that for 
the quadratic path-length dependence. 
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